Injection statistics simulator for dynamic analysis of noise in mesoscopic devices 
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We present a model for electron injection from thermal reservoirs which is applied to particle sim- 
ulations of one-dimensional mesoscopic conductors. The statistics of injected carriers is correctly 
described from nondegenerate to completely degenerate conditions. The model is validated by com- 
paring Monte Carlo simulations with existing analytical results for the case of ballistic conductors. 
An excellent agreement is found for average and noise characteristics, in particular, the fundamental 
unities of electrical and thermal conductances are exactly reproduced. 
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The systematic trend of reduction in the size of elec- 
tronic devices has led to the appearance of new phenom- 
ena that require special attention to be properly investi- 
gated. In particular, mesoscopic conductors are attract- 
ing increasing interest in recent years.ul3 Here, the mi- 
croscopic interpretation of carrier transport and fluctu- 
ations demands for approaches which differ from those 
typically used in macroscopic devices. Several techniques 
have been used to this end. Accordingly, to account for 
phase coherence the scattering matrix theory originallit 
proposed by LandauerH has been further elaborated.oEl 
A Wigner function formalism has also been uaad-for the 
analysis of ballistic and diffusive conductors.Q^Ei] When 
phase coherence does not play an essential role, semiclas- 
sical methods based on the Boltzmann-Lansavift equa- 
tion have shown to provide viable solutions .IlilEJ In all 
these theoretical approaches, the modeling of the con- 
tacts is proven to be crucial. The active region of the 
devices is considered to be surrounded by leads which 
are usually treated as ideal thermal reservoirs. In other 
words, contacts are assumed to be always at thermal 
equilibrium; absorbed carriers are thermalized immedi- 
ately, thus any correlation is destroyed, and emitted car- 
riers obey a Fermi-Dirac distribution. 

Recently, particle simulations, mainly based on the 
Monte Carlo method, have beerL— iised to study fluc- 
tuations in mesoscopic structures. lij^llS This technique, 
which is wides prea d for the analysis of macroscopic elec- 
tronic devices,ll3E3 has the advantage of being applicable 
under physical conditions which can be very far from 
thermal equilibrium, and often can not be studied ana- 
lytically. Moreover, it can provide detailed microscopic 
information about the physical processes and the time 
scales associated with transport and fluctuations in elec- 
tronic devices. This last feature makes particle simula- 
tions quite attractive for the study of mesoscopic con- 
ductors. In the case of macroscopic devices, the pres- 
ence of energy dissipation and diffusive regions inside 



the structures washes out the influence of the contact 
injecting statistics on the output currents and voltages, 
and the simulatipa of contacts does not require very de- 
tailed models .Ellu On the contrary, when dealing with 
mesoscopic structures, the modeling of carrier injection 
from thermal reservoirs is a delicate problem. In partic- 
ular, the statistics of electron injection associated with 
a Fermi-Dirac distribution at the electrodes is essential 
for the correct analysis of fluctuations and effects related 
to Pauli exclusion principle. For a classical iniector, the 
statistics of transmitted-jcharge is Poissonian,! and can 
be easily accounted for.Ej On the contrary, under degep^. 
erate conditions one should use a binomial distribution,^^ 
and to our knowledge this issue has not been addressed 
so far. 

The aim of this paper is to present a model for par- 
ticle injection from ideal thermal reservoirs into one- 
dimensional mesoscopic conductors which takes into ac- 
count the fluctuating occupancy of the incoming elec- 
tron states associated with a Fermi-Dirac distribution. 
The model can be continuously applied from nondegen- 
erate to degenerate statistics at the contacts. The re- 
sults obtained with a Monte Carlo simulation of a one- 
dimensional two-terminal ballistic conductor implement- 
ing the present contact model are compared with existing 
analytical results to validate the injection scheme. 

Let us consider a one-dimensional conductor connected 
to leads which act as perfect thermal reservoirs. The den- 
sity (in fc-space) of incoming electron states with wave 
vector k impinging per unit time upon the boundary be- 
tween the leads and the conductor, C^, is given by the 
product of the density of states n^ and the velocity Vk 
normal to the boundary, C,k — n^Vk — — — , where we 
have taken a parabolic isotropic e — fc relation. These 
rik states obey Fermi-Dirac statistics, thus only a frac- 
tion /(£fe) = {l-|-exp[(efe — ei?)//cBT]}~^ of them will be 
occupied and eventually will inject a carrier into the con- 
ductor, with ef the Fermi level. Therefore, the injection 



rate density of carriers with momentum k, Tk, is given 
by Ffc = Cfe/(£fe)- While (k does not depend on time, the 
instantaneous occupancy of an incoming fc-state /(sfc, t), 
of which /(efc) is the avpjage, fluctuates in time obeying 
a binomial distributionE3 with a probability of success 
/(efe). The injecting statistics imposed by this binomial 
distribution is determined by the Fermi-Dirac statistics 
electrons obey, i.e., ultimately, by Pauli principle. When 
Sk — £f ^ —ksT, f{ek) — 1 and the injecting statistics 
of the corresponding fc-state becomes uniform in time. 
On the contrary, when Sk — ep ^ ksT, f{sk) ^ 1 and 
the injecting statistics of the corresponding fc-state be- 
comes Poissonian in time. In a completely degenerate 
(quantum) reservoir, the former condition is fulfilled for 
any incoming fc-state and the injection is uniform in time 
for all the k values up to the Fermi wave vector kp. On 
the contrary, in a nondegenerate (classical) reservoir, the 
latter condition applies for all k values. 

To reproduce the injecting statistics imposed by the 
Pauli principle in a particle simulation, it is necessary 
to discretize momentum space into a certain number of 
meshes of width Afc around discrete values of fc, ki. For 
each of these meshes, the number of incoming electron 
states per unit time with wave vector ki is given by 
Cki,Ak = Cki^k, with a probability of occupancy given 
by /(cfci)- In the simulation, at each time interval of 
duration l/Cki,Ak an attempt to introduce an incoming 
electron with wave vector ki takes place. At this point a 
random number r uniformly distributed between and 
1 is generated, and the attempt is considered successful 
only if r < /(£&;) . This rejection-technique scheme prop- 
erly accounts for the injection statistics at each mesh in 
fc-space. 

For a completely degenerate reservoir, in every mesh 
up to fcp-| an electron is injected every time interval 
l/Cfei,AfcJ^ ^ii'i there is no need of the rejection tech- 
nique. This is the case of the simple contact modeling 
used in Ref. E4. For a nondegenerate reservoir, since 
fiski) ^ 1 for all fcj-states, it is possible to use a global 
Poissonian statistics characterized by an injection rate 
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/(J Tkdk. Accordingly, the time between two 
consecutive electron injections is generated with a prob- 
ability per unit time P{t) = FcZase~^'''"=*. Then, the 
electron wave- vector is randomly picked from a Maxwell- 
Boltzmann distribution, and there is no need of using a 
mesh in fc-space. This is the scheme used in Refs. 15-1^. 



For any intermediate level of degeneracy, to account for 
the proper statistics at each value of ki, it is necessary to 
use the scheme explained above, which of course is also 
valid in the classical and degenerate limits, but less effi- 
cient from the point of view of computation time. The 
accuracy of the proposed scheme depends on the number 
of meshes in fc-space used to inject carriers. In any case, 
it is not necessary to use a very large number of meshes. 
Indeed, it is well known that the noise of an electrical 
system depends only on the kinetics of electron states in 
a small energy range around the Fermi level. We have 
checked that very satisfactory noise results can be ob- 
tained by simulating just an energy range of Sfc^T above 
and below ep, and dividing it into 50 meshes. Below and 



above this range all fc-states can be respectively consid- 
ered to be completely occupied and empty. Therefore, 
they do not contribute to current fluctuations and can 
be ignored. 

In the following we will report the results of a Monte 
Carlo simulation of a one-dimensional two-terminal bal- 
listic conductor of length L connected with two thermal 
reservoirs modeled according to the above scheme. The 
temperature is taken to be 300 K and the effective mass 
m = 0.25 Too, toq being the free electron mass. Carriers 
are considered to move ballistically into the conductor 
following the classical equations of motion, and when a 
voltage U is applied to the leads, electrons are accel- 
erated by an electric field E = U/L. When a carrier 
inside the conductor reaches a contact, it is considered 
to be immediately thermalized and it is cancelled from 
the simulation. For simplicity, the cross-sectional area 
of the conductor is assumed to be sufficiently small so 
that only the lowest sub-band is occupied. We remark 
that Coulomb interactions are ignored and Pauli exclu- 
sion principle is taken into account only at the contact 
injection. Under these conditions, the literature provides 
several analytical results which will be used to validate 
the model. 

Figure H shows the low-frequency value of the current 
spectral density Sj (0) at equilibrium normalized to 2qls 
as a function of the degeneracy factor ep/ksT, with sp 
measured with respect to the bottom of the conduction 
band. Is = Q Jq Vknkf{£k)dk is the saturation current, 
i.e., the maximum current a contact can provide. In the 
classical limit, corresponding to large negative values of 
the degeneracy factor, S'/(0) — Aqls- Here all carriers 
contribute to the current noise and Si{0) is just the sum 
of the full shot noise related to th&Jsmx opposing cur- 
rents Is injected by the contacts.QBESEa Under degen- 
erate conditions, corresponding to positive values of the 
degeneracy factor, Si{0) decreases with respect to iqls 
in accordance with the suppression factor e^ksT re- 
lated to Fermi correlations at the reservoirs JjiO Here, as 
known, only carriers around the Fermi level contribute 
to the noise. As shown by the figure, the agreement be- 
tween the results of the Monte Carlo simulation and the 
analytical expectations in the nondegenerate and degen- 
erate limits is excellent, thus indicating that the carrier 
injecting statistics achieved with the proposed model is 
valid in both regimes. 

Figure reports the current autocorrelation function 
Ci{t) in the same structure of Fig. 1 when ep/ksT = 100 
for several applied voltages as a function of time nor- 
malized to the transit time at the Fermi level tt — 
L{ra/2epY^'^ . The corresponding I — U curve is shown 
in the upper inset of Fig. ^ Due to the unbalance in 
the number of carriers reaching the opposite contact, the 
current increases linearly with the applied voltage un- 
til qU = ep, when the current reaches the saturation 
value Is = 2q£p/h. For higher U all the carriers in- 
jected at the left lead reach the opposite contact, while 
no electron injected at the anode reaches the cathode, 
and therefore the current saturates. The conductance in 
the linear region corresponds to the value of the funda- 
mental unit 2q^ /hn The current autocorrelation function 



C/(i) exhibits the following features. At equilibrium it 
shows the typical triangular shape, vanishing at t = tt 
since the contributions of carriers moving in both direc- 
tions are symmetrical. This shape parallels that of a 
vacuum tube with a constant velocity emitter. Here, the 
same shape comes from Pauli principle which allows only 
carriers in a small range around the Fermi energy, and 
therefore moving with practically the same velocity, to 
contribute to the noise. When a voltage below ep/q is 
applied to the structure, Ci{t) exhibits a two slope be- 
havior because now the transit times of carriers moving 
in opposite directions are different. At voltages higher 
than ep/q a negative part appears in C/(i) since the car- 
riers injected against the electric field no longer reach 
the cathode and return to the anode. At further increas- 
ing voltages the negative part appears sooner due to the 
shorter time it takes to the carriers to return back to the 
right contact .0 The low- frequency spectral density 5/(0) 
as a function of U is shown in the lower inset of the same 
figure. At equilibrium the value obtained corresponds 
to 5*7(0) = Sq^ksT/h, which, when compared with the 
Nyquist formula 5/(0) = AksTG, provides again for the 
static conductance G — 2q^ /h. The fact that in our 
model the conductance obtained from the I — U curve 
reproduces the fundamental unit value is a valid check of 
the correct use of the one-dimensional density of states 
at the contacts. Furthermore, the same value 2q^/h is 
also obtained from noise results at equilibrium, which 
proves that the injecting statistics model under degener- 
ate conditions here proposed is also correct. The voltage 
dependence of 5/(0) exhibits a step- like behavior, taking 
the equilibrium value up to U — sp/q and half this value 
for higher U. This behavior is understood as follows. 
For U < Sp/q, carriers around the Fermi level which 
are injected at both contacts reach the opposite side and 
therefore both contribute to the low-frequency noise. On 
the contrary, for U > sp/q only carriers injected at the 
cathode reach the anode and thus the value of 5/(0) is 
halved. We remark that all the results shown in Fig. g are 
in excellent agreemeal-with previous analytical results in 
degenerate systemsE a 

One of the advantages of using particle simulations for 
the noise analysis is the possibility to interpret the time 
and frequency behavior of fluctuations in terms of dif- 
ferent contributions. Thus, in Fig. ^, in the case of de- 
generate conditions and qU/ep — 1.01, C/(i) is decom- 
posed into velocity Cy(t), number CN{t), and velocity- 
number CvN{t) contributions.Ej Here it can be observed 
that the origin of the negative part in C/ (t) comes from 
the velocity-number correlation, which at zero time ex- 
actly compensates the velocity contribution. Further- 
more, from the time dependence of fluctuations three dif- 
ferent characteristic times can be identified. The shortest 
one corresponds to the transit of carriers injected from 
the left contact, as better evidenced in Cn and Cvn- It 
is close to the transit time at equilibrium tt, but slightly 
shorter due to the acceleration of the field. A second one 
is the time taken by the carrier injected at the anode to 
reverse its velocity, and is reflected mainly in Cvn- Fi- 
nally, the longest one is the time at which all correlation 
functions vanish. It corresponds to the time spent by 



the electrons injected at the anode to return back to the 
same contact. 

As final result, in Fig. |j we report the spectrum of 
the thermal conductivity k(/) at equilibrium, calculated 
according to Ref. nfl. Here, we remark that not only the 
correlations of electrical current fluctuations are involved, 
but also those of the heat flux and the cross-correlations 
between both. The oscillatory structure of Re[K(/)] , with 
geometrical resonances at the inverse of tt, is associated 
with the fact that all the involved correlation functions 
exhibit the triangular shape already found in the case of 
the current (see Fig. g). A corresponding structure is 
detectable also in the imaginary part Im[K(/)]. Again, 
the results of the Monte Carlo simulation are in excellent 
agreement with analytical results. In particular, they 
reproduce with great accuracy the fundamental unit of 
thermal conductance K — 2Tr^k^T/3h inferred in Ref. 
\IX\ , which again conflrms the validity of the proposed 
model for the injection statistics at the reservoirs. 

In summary, we have presented an injection scheme of 
electrons at thermal reservoirs for particle simulations of 
mesoscopic conductors which takes into account the bi- 
nomial distribution of the injected electrons imposed by 
Fermi statistics. The model has been validated for the 
case of ballistic transport in quasi one-dimensional degen- 
erate conductors. In particular, at thermal equilibrium 
we have reproduced the fundamental units of electrical 
and thermal conductances, analytically calculated from 
the correlation-function formalism. The scheme can be 
applied continuously from classical to completely degen- 
erate conditions, and it can be extended to two and three 
dimensions, and multi-subband systems. The proposed 
scheme is open to applications involving degenerate dif- 
fusive conductors. 
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FIG. 1. Low-frequency current spectral density at equilib- 
rium normalized to 2qls as a function of the degeneracy factor 
ef /ksT, at T = 300 K. Symbols correspond to Monte Carlo 
calculations and lines to nondegenerate (solid) and degenerate 
(dotted) conditions, respectively. 



FIG. 2. Current autocorrelation function normalized to the 
zero-time value for several applied voltages. The calculations 
correspond to degenerate conditions with ep/kBT — 100. 
The upper inset reports the I — U characteristic of the struc- 
ture, while the lower inset shows the 5/(0) — U characteristic. 
Results are normalized to appropriate units. Solid lines in 
the insets correspond to analytical results and full circles to 
Monte Carlo calculations. 



FIG. 3. Decomposition of the current autocorrelation func- 
tion into velocity, number and velocity-number contributions 
for the case of qU /ef ~ 1.01 and ef /ksT = 100. 



FIG. 4. Spectrum of thermal conductivity at equilib- 
rium normalized to the theoretically predicted zero- frequency 
value. Solid and dotted lines correspond to the real and 
imaginary parts, respectively. Calculations are performed for 
EF/ksT = 100. 
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